UG-FT-82/97 
November 1997 



Supersymmetric calculations with component fields 
in differential renormalization a 



F. del Aguila, 1 A. Culatti, 1,2 R. Munoz Tapia, 1 and M. Perez- Victoria 1 
1 Dpto. de Fisica Tedrica y del Cosmos, Universidad de Granada, 
18071 Granada, Spain; 
Dip. di Fisica, Universitd di Padova, 35131 Padova, Italy 

Two applications of the method of differential renormalization to superymmetric 
gauge theories are reviewed. The photon propagator in supersymmetric QED is 
renormalized at one loop and the first supergravity contributions to the anomalous 
magnetic moment of a charged lepton are obtained. 



1 Introduction 

Quantum Field Theories are in general plagued with infinities, which are usu- 
ally regularized with some sort of cut-off. Then the renormalization program 
allows to get rid of the singularities in a consistent way by a redefinition of the 
parameters in the Lagrangian. In some methods, like BPHZ, renormalization is 
carried out in one single step without intermediate regularizationEI. Whatever 
scheme is used, it is desirable that it maintains the symmetries of the theory 
under study. Theories with a high degree of symmetry can be very demanding 
in this respect. This is the case of gauge and supersymmetric theories. In a 
scalar theory with very little symmetry, almost any regulator does the job. 

Dimensional regularizationEI is a remarkable method that explicitly pre- 
serves gauge invariance. Because of its relative simplicity and the importance 
of gauge theories in the description of Nature, this method has become stan- 
dard in many computations in Quantum Field Theory. However, dimensional 
regularization has difficulties in coping with supersymmetric theories: it can 
break supersymmetry (SUSY). Essentially, the reason is that the equality of 
Bose and Fermi degrees of freedom only holds for specific values of the space- 
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time dimension, which is changed by this regularization. There exists a variant 
of the method, dimensional reduction 13, in which the field components remain 
unchanged, while the loop integrals are performed in a d-dimensional space. 
This method preserves SUSY, at least at lower orders, but at higher ones the 
situation is more involved, especially in the case of broken SUSY (see Ref. □ 
and references therein) . Another approach is based on higher derivative reg- 
ularization supplemented by Pauli-Villars What is clear in any case is the 
scarcity of simple methods suitable for supersymmetric theories. 

The method of differential renormalization (DR) has appeared rece ntlyi. 
It is a method of renormalization without regulators or explicit counterterms 
which works in coordinate space. DR does not modify the space-time dimen- 
sion, thus being a candidate for preserving SUSY. The method has been shown 
to be quite powerful in a three-loop computation for the massless A</> theory. 
Other applications include lawer dimensional theories, Chern-Simons, naoper- 
turbative calculations, etc. 0. Formal aspects like checks of unitaoty 119, its 
relation with dimensional regularization J_hc inclusion of masses u and the 
consistency of the procedure to any order liil have-been also studied. Differ- 
ent versions of DR have been developed in Ref. The first application to 
supersymmetric theories-was a calculation of the /3-function to three loops in 
the Wess-Zumino modelE^I. In Ref.LLa DR was applied to pure supersymmetric 
gauge theories, and the /3-function was obtained to two loops and one loop, in 
the abelian and non-abelian case, respectively. These calculations were per- 
formed with superspace techniques. Very recently this method has also been 
employed in non-perturbative calculations in supersymmetric gauge theoriesuJ. 
In the superspace formalism, SUSY is manifestly preserved; however the situa- 
tion is more involved in the physically interesting case of broken SUSY, where 
it is usually preferred to work with component fields (see, however, Ref. Ej). 
Here we review the two existing applications of DR to supersymmetric gauge 
theories in the component approach: the vacuum polarization in supersymmet- 
ric QED (SQED)liZl and the calculation of the anomalpns magnetic moment of 
a charged lepton, (g — 2);, in unbroken supergravityEj. These two examples 
provide non-trivial tests of the potential of DR to preserve (abelian) gauge 
invariance and SUSY. In standard DR, the fulfilment of the corresponding 
Ward identities is accomplished adjusting at the end the different scales that 
appear in the renormalization procedure. In Ref. E3 we described a procedure 
to constrain the scales appearing in DR in such a way that Ward identities are 
automatically satisfied. We_have verified the consistency of this approach in 
different one-loop examples tL3. Here we shall use this constrained procedure of 
DR. 

In what follows we briefly review the method of DR and its constrained 
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form. Then we present the two one-loop supersymmctric calculations with 
component fields. The first example, the renormalization of the photon prop- 
agator in SQED, illustrates the method. Afterwards the more involved case of 
the evaluation of (g — 2)i in unbroken supergravity is revised, and the results 
are discussed. The last Section is devoted to conclusions. 

2 Differential renormalization and its constrained version 

In coordinate space the amplitudes are finite as long as point coordinates are 
kept apart; singularities only arise when points coincide. These singularities, 
when too severe, give rise to divergences in integrals on internal points or in 
Fourier transforms. In other words, the bare expressions are in general ill- 
defined distributions. The idea of differential renormalization is to rewrite 
these expressions as derivatives of less singular ones. The derivatives are un- 
derstood in the sense of distribution theory, i.e., acting formally by parts. The 
amplitudes written in this way are identical to the bare ones for separate points 
but behave well enough at coincident points. One has then the renormalized 
amplitudes in coordinate space. In this process some dimensionful integration 
constants appear that will play the role of renormalization scales. An illustra- 
tive example is given in massless Xcj) 4 theory, where the four point function at 
one loop has a factor 1 /x 4 (x being the difference of any two external points) . 
The renormalization is given by the substitution 

1_ I log(M 2 x 2 ) 

x* 4 x 2 [ ' 

The scale M is the renormalization scale alluded to above. Notice that the 
r.h.s. in Eq. ([!]) has a well defined Fourier transform. We will work in Euclidean 
space, where the handling of the expressions is.simpler. Analogous equations 
can be obtained for other singularities (see Ref.O for details) and the_program 
can be in general carried out for any theory at any number of loops til. 

In principle, the scales introduced for different diagrams are independent. 
These scales can be choosen in different ways (as long as the same scale appears 
in identical subgraphs), each corresponding to a choice of the renormalization 
scheme. The easiest one is to take all scales to be equal. This choice is ap- 
propriate for a scalar \<p 4 theory, for instance. However, in theories with more 
symmetry much care must be taken when fixing the constants. In gauge theo- 
ries for example, there are restrictions to be fulfilled which are dictated by the 
Ward identities. The usual way of proceeding is to renormalize the amplitudes 
and then relate the constants imposing the Ward identities. Although this is a 
possible procedure, one would prefer that gauge symmetry were automatically 
preserved, as occurs in dimensional regularization. 
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The constrained version of DR fixes the arbitrary constants from the be- 
ginning (see Ref. Ell for details). The idea is to restrict oneself to a set of 
consistent rules to manipulate the singular expressions. Besides the usual DR 
rules, differential equalities and formal integration by parts, two other condi- 
tions appear to be enough to conveniently fix the local terms at one loop. First, 
the factorization of delta-functions in the renormalization procedure, i.e., 

[F(x, xx,..., x n )6(x - y)] R = [F(x, x x ,..., x n )] R 8{x - y), (2) 

where F is an arbitrary function and R stands for "renormalized" . Second, 
one demands that the propagator equation (□ — m 2 )A m (a;) = — 5{x) has a 
general validity when embedded in any amplitude. This is to say that 

F(x,xi, . . .,£„)(□ - m 2 )A m (x) = F{x,x\, . . . ,x n )[-8(x)] (3) 

holds for F arbitrary, and not only for well-behaved enough functions. Above, 
A m (x) = 4^2 mKl ^ mx ^ is the Feynman propagator for a particle of mass m, 
with K\ a modified Bessel function. These rules allow to renormalize a set of 
basic functions which are used to expand the amplitudes. The basic functions 
for tadpole, bubble and triangular diagrams with massless propagators are 
denoted by 

A = A(x)5{x) , (4) 
B[0] = A(x)O x A(x) , (5) 
T[0} = A(x)A(y)G*A(x-y) , (6) 

where A(x) = -^i^i is the massless propagator and O x is a differential oper- 
ator. For example, using these rules one can easily obtain 

T R [a] = -B R [l}(x)6(x-y) 

1 1 \ogx 2 M 2 „ 

The massive counterparts are defined accordingly and will be denoted with a 
subindex m. The DR substitutions for massive expressions canj-he obtained 
using recurrence relations among Bessel functions (see Refs.El andli3 for a more 
detailed discussion). The only one we shall need is 



m 2 K 1 (mx) 2 1 2 mK (mx)K 1 (mx) 2 M 2 

2 * ~Z y— ' 4?71 ) -|- 7T log g- 

x^ 2 x m 



where M = =^-, and 7r = 1.781... is Euler's constant. When one has to handle 



tensor structures, the general procedure involves decomposing the expressions 
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into trace and traceless parts. Consistency with Eqs. (|^, |^) may imply the 
appearance of local terms. For instance, in the renormalization of T[<9 M <9„] one 
has 

T R [d,d v ] = i VT*[D] + T[d„d„ - ±8^n] -±^ S (x)6(y)6 tlu . (9) 

The last term is the local term alluded to, which is needed for the consistency 
of the equality: 

B R [d,](x)S(y) = -^□»T[1] + □"TPJ - 2fl^T[0 <r ] d*T R p] 

+ 2d v a T R [d^]+T R [d^D] . (10) 

To complete the renormalization of T[d M <9„], one must note that only the first 
term in Eq. (^) is singular and requires renormalization (see Eq. (0)). Eq. (^|) 
shows that in general 

<WT R [9,A] ^ [<WT[fl M $,]] fl . (11) 



3 Vacuum polarization in SQED 

This is a simple example of a calculation in a supersymmetric abelian gauge 
theory using component fields. As a matter of fact we will only check the 
transversality of the photon self-energy and not any SUSY relation. The vac- 
uum polarization in SQED has two contributions, one coming from the scalar 
loops (the corresponding diagrams are depicted in Fig. 1), and the other from 
QED. Before renormalizing, we will reduce their expressions to sums of basic 
functions and add. the two parts. The Feynman rules in coordinate space can 
be found in Ref.lia. 

The vacuum polarization in scalar QED reads 

n$(a:) = -e 2 A m {x) d^dv A m (x) , (12) 
11$ (x) = -2e 2 6„ v A m (x)6(x) , (13) 

where A d B = A9B — B9A. In terms of basic functions one has 

n£)(a:) = -e 2 {4B m [9 Ai a,] - d^B m [l}} , (14) 

11$ (x) = -2e%„A m . (15) 

Now the propagator equation can be used to rewrite the tadpole contribution, 
Eq. (|l5|), in terms of bubble functions, 

A m = A m (x)5(x) 
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Figure 1: One-loop diagrams contributing to the vacuum polarization in scalar QED. 

= -A m [x)(a -m 2 )A m (x) 

= -B m [D]+m 2 B m [l] . (16) 

On the other hand the QED contribution (diagram 1 in Fig. 1, but with a 
fcrmion running in the loop) can be written 

IVOz) =4e 2 {(m 2 V + ^^D-9 M ^)B m [l]+2B TO [^]-VBmp]}- (17) 

Then the supersymmetric QED vacuum polarization is the sum of the spinor 
diagram plus twice the scalar ones. One directly obtains a transverse result 
depending only on one basic function: 

Tl^(x) = -2e 2 {d^ - ^□)B ro [l] . (18) 

This equation exhibits the consistency of constrained DR and abelian gauge 
invariance in this simple one-loop supersymmetric calculation. The renormal- 
ization is completed substituting B m [l] by its renormalized expression, given 
by Eq. (§: 

__d , , e 2 ,„ . ... mKn(mx)Ki(mx) 

u U x ) = -j^fi 9 ^ - <V D )[( n - 4™ 2 ) — ov ; u — L 

M 2 

+ 2^ 2 log — 6(x)} . (19) 
m 



6 



4 Supergravity contributions to (g — 2); 



This calculation is more involved. In order to show how coordinate space tech- 
niques can be used to calculate a quantity typically defined for fixed momenta, 
we first work out the standard QED correction in coordinate space. 

4-1 (g — 2)i in momentum and coordinate space 

The anomalous magnetic moment of the electron, (g — 2)/, is usually defined in 
momentum space as a static limit. For p and p' being the incoming momenta of 
the electrons and q = p + p' being the outgoing momentum of the photon, the 
parity conserving vertex containing all radiative corrections can be expressed, 
for on-shell external electrons, in terms of two form factors: 

A^pi) = t e{(F 1 ( q 2 ) + F 2 ( q 2 )h, + ^^(p,-p>J} . (20) 
Then the anomalous magnetic moment is defined as 




The corresponding expression for the vertex in coordinate space is related to 
the momentum space one via 

A M (p,p') = J d A xd 4 ye^ x e l P'-yA^x,y) . (22) 

The on-shell condition for the external electrons corresponds to imposing Dirac 
equation on their wave functions. This means that we must substitute terms 

of the form ^ x f{x 1 y) and f(x,y) y in A fl (x,y) by —mf{x,y) and mf(x,y), 
respectively. Similarly, the static limit corresponds in coordinate space to 
imposing UA^ = 0. Therefore, terms in A fl (x,y) like O z f(x,y) = [—(d* + 
®%,)][~ (^u + @%,)]f( x >y) must be neglected. In practice, however, one cannot 
always extract all derivatives, and pieces containing internal derivatives remain. 
Hence, Dirac and Maxwell equations cannot be directly imposed. The simplest 
solution is to perform at the end a Fourier transform to obtain A^(p,p') and 
to take the appropriate limits. This procedure may seem to end up in the 
usual momentum space one. The difference, however, is that renormalization 
is carried out before Fourier transforming, so the Fourier integrals are finite 
and can be computed in four dimensions without any regulator. 

In QED, (g — 2)i vanishes at tree level (this is predicted by Dirac's equation: 
g = 2 for a particle of spin 1/2), but it is well-known that a finite non-zero 
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(g — 2)i is generated at one loop. At this order the vertex correction, V^, is 
given by the standard triangular diagram and reads 

Vpfay) = He) 3 7a (r -m)A m (x)^(- 0» - m)A m {y) la A{x - y). (23) 



Using Leibnitz rule to rearrange derivatives, Eq. ([23|) can be expressed in terms 
of the triangular functions defined in Eq. ([|), in this case with two massive 
propagators, i.e., 

T m [0] = A m (x)A m (y)O x A(x - y) . (24) 

Using Dirac's equation whenever possible the resulting expression is 

V^x, y) = *e 3 {-4 lfl d* ■ 0»T m [l] + 4(9* - d«)T m [0\ - AjM ~ d v a )T m [d a ] 
+ 4mT m [a M ] + 27 M T fn [D] - *y o T m [0 o M ]}. (25) 

The terms proportional to 7p do not contribute to {g — 2); and we ignore them 
in the following. Of the rest, only the last triangular function, T m [<9 a <9 M ], is 
singular. As in Eq. (^J) it is decomposed into a part proportional to 5^ v and 
a finite traceless part. Only this traceless part contributes to (g — 2);. Hence, 
(g — 2)i can be extracted from the finite expression, 

V M (x, y) g - 2 = *e 3 {-4(^ - 0£)T m [fl + 4mT m [d M ] - 4 7a T m [d Q ^ - j^D]} , 

(26) 

which can be readily Fourier transformed (see Eq. (f22[)). The necessary inte- 
grals in the static limit, 



are 



V A = V = -m 2 , <f , (27) 



T m [^] = -^4^2^-K), (28) 



327r 2 m 2 

1 i 1 1 

T m [9 Q 9 M - -rS afl D] = - 2 {-^{PaP'+P'aP^) + ~(PaPn + p' a P 



+7™ 2 <W} , (29) 



Now the value of the anomalous magnetic moment can be read from the coef- 
ficient of p M — p'^ in U M (a;,y) 9 ~ 2 (see Eq. (po|)): 

<7 — 2 2m ie 3 1 L a , . 

2 ze loTT^m 2 4 27r 

This is the well-known Schwinger result 111. 



Dl D2 D3 D4 D5 




D6 D7 D8 D9 DIO 



Figure 2: Diagrams of order en 2 contributing to (g — 2); in supergravity. A graviton is 
exchanged in diagrams D1-D5 and a gravitino in D6-D10. 

4-2 (g — 2)i in unbroken supergravity 

In a supersymmetric theory (g — 2)/ vanishesJiecause no such term appears 
in the Lagrangian of a chiral supermultiplet EJ. Hence, as long as SUSY is 
preserved, all quantum corrections must cancel order by order. Therefore, the 
anomalous magnetic moment of the lepton besides being an observable, is also 
an ideal arena to check theoretical implications and to perform consistency 
tests of regularization methods. Ferrara and Remiddi also proved explicitly 
that in global SUSY the one-loop QED corrections, order e 3 , do cancel. This 
is to say that the fermion contribution (Schwinger result) and the scalar one 
cancel each other. The latter results from twice the same triangular diagram 
but with the slepton and the photino replacing the lepton and the photon, 
respectively. 

The one-loop gravitational corrections_are of order en 2 — %ireGN, result- 
ing from a graviton or gravitino exchange Ej (the corresponding diagrams are 
depicted in Fig. 2). Using dimensional regularization, Berends and Gast- 
mansE3 calculated the five diagrams where a graviton is exchanged. All five 
diagrams are infinite but their sum is finite. The finiteness of (g — 2); in a 
nomjenormalizable theonz_such as gravity seemed miraculous. Del Aguila et 
al. t3 and Bellucci et al. Ej checked that when gravitation is embedded in a 
supersymmetric theory (unbroken), the contributions from_the graviton and 
the gravitino cancel, as required by SUSY. Bellucci et al. also traced back 
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to an effective chiral symmetry in the gravitino sector the finiteness of the 
gravitino contribution and then of the graviton contribution, if their sum has 
to vanish. Dimensional regularization does not yield a vanishing value for 
(g — 2) i. This is one example of a case where dimensional regularization breaks 
SUSY. A (one-loop) SUSY preserving method is required in order to obtain 
such a cancellation and this was shown to be the case in dimensional reduc- 
tion. In Ref. Ea we calculated these contributions using DR (the Lagrangian, 
Feynman rules and other technical details can be found there). The use of 
the rule extending the validity of the propagator equation allowed us to relate 
diagrams with different topology before explicit renormalization. Then consid- 
ering the sum of each diagram and its supersymmetric partner, we found that 
the contibution of each sum to (g — 2); depends only on one singular scalar 
basic function, T m [d], apart from finite terms. (This is analogous to what 
happens in the calculation of the vacuum polarization in SQED above.) The 
terms proportional to T m [D] cancel out in the complete sum and so do the 
finite terms. In this way SUSY is preserved, i.e., a vanishing value of (g — 2); 
is obtained. Terms proportional to did not appear, either. Hence, on-shell 
gauge invariance, forbidding these terms, is also respected in this calculation. 
It is worth to note, however, that singular basic functions with other tensor 
structure do appear in individual graphs and in the total graviton (gravitino) 
contribution. Then, for these contributions which become physically meaning- 
ful if SUSY is broken, the tensor decomposition of these singular functions and 
the local terms they introduce (see for instance Eq. (^) are relevant. In Ref.tLa 
we used the engineering trace-traceless decomposition neglecting local terms. 
We call this method 'partially constrained' DR. Therefore, the partial results 
may and do differ when using 'partially constrained' DR or constrained DR 
(for the latter includes the local terms), although the total sum is the same. In 
Table ^ we gather the results in dimensional regularization, dimensional reduc- 
tion, 'partially constrained' DR and constrained DR. As can be observed, the 
last three methods (columns) preserve SUSY (add to zero) , but 'partially con- 
strained' DR gives a different value for the graviton (gravitino) contribution. 
The graviton (gravitino) contribution depends on the regularization method 
and is not well defined by itself. This is related to the fact that we are deal- 
ing with a non-renormalizable theory, and in the non-supersymmetric case no 
symmetry protects the value of {g — 2);. The fact that dimensional reduction 
and constrained DR give the same result for this contribution seems to indi- 
cate that the requirement that the regularization/renormalization method be 
compatible with gauge invariance and SUSY greatly constrains the results at 
one loop. 
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Table 1: Contributions of the diagrams in Fig. 2 to ( g 2 2 ) in units of Gn ™ , obtained 
with dimensional regularization, dimensional reduction, 'partially constrained' DR and con- 
strained DR. 



Diagram 



Dimensional 
Regularization 



Dimensional 
Reduction 



'Partially 
Constrained' DR 



Constrained DR 



Dl 

D2+D3 

D4+D5 

Graviton 

(D1+D2+D3 

+D4+D5) 



1 1 _ 61 
3 n-4 36 



1 1 _ 29 
3 n-4 18 



11 1 
3 n-4 



32 
!) 



11 1 

3 n-4 



35 
!) 



+ 7 

n — 4 



7/4 



-4—4 + 6 

n — 4 



.I log f Mil - 25 

6 10S V m2 y 18 



1/2 



1/2 



D6 



D7+D8 



D9+D10 

Gravitino 

(D6+D7+D8 

+D9+D10) 

TOTAL 

(Graviton 

+Gravitino) 



8 1 _ 55 
3 n-4 18 



4 1 _ 13 
3 n-4 9 



-4-4+4 

n — 4 



8^L 37 

3 n-4 18 



4 1 _ 4 
3 n-4 9 



-4^-r +2 

n — 4 



log 



( Ml) + la 

y^r j ^ is 



■ log 



■ log 



(£)+« 



- 1 



-1/2 



5/4 



-1/2 



-1/2 



5 Conclusions 

DR is a method of renormalization recently proposed!, which works in coor- 
dinate space and does not introduce any intermediate regulator. It seems to 
have the potential of preserving gauge and chiral invariance. This procedure 
has been applied in several contexts and in particular to several supersym- 
metric calculations with component fields, which we have reviewed here. The 
method can be constrained to get rid of arbitrary constants, except for the 
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neccesary renormalization group scale. Constrained DR is defined by a set of 
rules which completely fix the renormalization of singular expressions at least 
to one loop. In the two examples worked out, a transverse renormalized vac- 
uum polarization in SQED and a vanishing (g — 2)i in supergravity have been 
obtained. The use of the rule extending the validity of the propagator equation 
plays an essential role: it allows to relate the expressions appearing in different 
graphs and enforce the supersymmetric and gauge invariance constraints from 
the beginning. In both cases, SUSY cancellations make the result insensitive 
to the renormalization of basic functions with non-trivial tensor structure, and 
thus to the inclusion of local terms in the tensor decomposition. From only 
these two calculations, however, it cannot be said how general this effect is. At 
any rate, if SUSY, is broken it is clear that the local terms become relevant. As 
shown in Ref. 113, the extensive use of the rules completely determines these 
terms, so constrained DR can in principle be used in the case of broken SUSY. 
Here we have only presented a simple (but non-trivial) consistency check of 
the constrained DR method for supersymmetric abelian gauge theories. A real 
test (or proof) should consider the Ward identities of both SUSY and gauge 
invariance, as well as higher orders. 
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